Bi-partite and global entanglement in a many-particle system with collective spin 

coupling 
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Bipartite and global entanglement are analyzed for the ground state of a system of N spin 1/2 par- 
ticles interacting via a collective spin-spin coupling described by the Lipkin-Meshkov-Glick (LMG) 
Hamiltonian. Under certain conditions which includes the special case of a super-symmetry, the 
ground state can be constructed analytically. In the case of an anti-ferromagnetic coupling and for 
an even number of particles this state undergoes a smooth crossover as a function of the contin- 
uous anisotropy parameter 7 from a separable (7 = 00) to a maximally entangled many-particle 
state (7 = 0). From the analytic expression for the ground state, bipartite and global entangle- 
ment are calculated. In the thermodynamic limit a discontinuous change of the scaling behavior 
of the bipartite entanglement is found at the isotropy point 7 = 0. For 7 = the entanglement 
grows logarithmically with the system size with no upper bound, for 7 7^ it saturates at a level 
only depending on 7. For finite systems with total spin J = N/2 the scaling behavior changes at 

7 = 7crit = l/J. 

PACS numbers: 03.65.Ud,03.67.Mn,73.43.Nq 
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I. INTRODUCTION 

Since the early days of quantum theory it was real- 
ized that quantum systems can possess correlations that 
do not have a classical counterpart 0, 0, El • F° r a l° n g 
time this phenomenon, called entanglement, has been of 
interest mostly in the context of foundations of quan- 
tum mechanics. With the advent of quantum information 
science it has been realized that entanglement is an 
essential resource for efficient computation and commu- 
nication. This initiated a more systematic study of its 
properties. While entanglement in small systems is by 
now well understood, many particle entanglement is still 
a widely open field. It is well known that quantum corre- 
lations and entanglement naturally occur in interacting 
many-particle systems but we only begin to understand 
their role in the systems 

Recently the entanglement properties of quantum sys- 
tems near the critical points of quantum phase transitions 
have attracted much attention. Two-particle entan- 
glement has been studied in terms of the concurrence H 
e.g. in one dimensional spin chains 
The concurrence contains however only limited informa- 
tion about the global distribution of entanglement, and 
other measures such as the bi-partite entanglement be- 
tween blocks of spins may be of larger interest. Bi-partite 
entanglement has been analyzed e.g. in one dimensional 
quantum spin chains in |12j . where it was shown that it 
scales logarithmically with the system size at the criti- 
cal point with a prefactor determined by the universality 
class and saturates in the noncritical regime. 

We here study the bi-partite entanglement in a system 
of spins with a collective coupling described by a gen- 
eralization of the Lipkin-Meshkov-Glick (LMG) Hamil- 
tonian [l3j |. Since this Hamiltonian is symmetric un- 
der the exchange of particles the Hilbert space sepa- 



rates into subspaces whose dimensions grow only linearly 
with the number of particles, which makes them numer- 
ically accessible. In particular we consider the case of 
an even number of spins and an anti-ferromagnetic cou- 
pling. Furthermore we concentrate on the case where 
the Hamiltonian can be factorized in a product of two 
terms each being linear in the total spin operators. Un- 
der these conditions, which includes the special case of a 
super-symmetry (SUSY) jl3, 0, the ground state can 
be constructed explicitly |16| . This ground state under- 
goes a smooth transition from a separable to a maximally 
entangled state as a function of a parameter 7, which 
characterizes the asymmetry between the collective spin 
coupling in x and y directions. The two-particle concur- 
rence in this system, which due to symmetry is the same 
for all pair of spins, has been analyzed in |l7| . Very re- 
cently also bi-partite entanglement has been studied in 
|18|| and |l9l | , however for the ferromagnetic version of the 
LMG model, where there is a quantum phase transition. 
Although in the anti-ferromagnetic case, considered here, 
there is no quantum phase transition, we find a discontin- 
uuous behavior of the entanglement when 7 is changed. 

After discussing the LMG model and its ground state 
under the condition of super-symmetry in Sec. II, we ana- 
lyze the bi-partite entanglement of this state in terms of 
the von Neumann entropy |2(j in Sec. III. We show that 
for an isotropic interaction (7 = 0) and in the thermody- 
namic limit the entropy grows logarithmically with the 
spin of the subsystem with no upper bound. On the other 
hand for any non vanishing 7 the entropy has an upper 
limit determined solely by 7. Furthermore it becomes a 
function of the ratio of the subsystem spin to the total 
spin rather than a function of the subsystem spin alone. 
In a finite system the transition between isotropic and 
anisotropic behavior occurs at 7 cr it = J . To under- 
stand the saturation of the entanglement quantitatively, 
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we give in Sec. IV an analytic estimate for the global 
entanglement by determining the geometric measure of 
entanglement. 

II. COLLECTIVE SPIN COUPLING AND 
SUPER-SYMMETRY 

Let us consider an even number N of spin 1/2 particles 
interacting through a nonlinear coupling of the collective 
spin J M = YljLi &°n> where <r M denotes the /i's component 
of the single-particle spin. The interaction is assumed to 
be of second order in the total spin and is thus a general- 
ization of the Lipkin-Meshkov-Glick (LMG) Hamiltonian 


H = aJ z + f3J 2 + J 2 - 2nJ v . (1) 

a and (3 are positive real numbers and thus the coupling 
is of the anti-ferromagnetic type. H commutes with the 
total spin J 2 and thus the total Hilbert space separates 
in sub spaces determined by the spin quantum number 
J. We here restrict ourselves to the case of maximum 
spin, i.e. J = N/2. As has been shown in [T^j, can be 
written as a product of two terms linear in the collective 
spin operators if j3 = a 2 : 

H = (aJ x + iJ y — i/ij (aJ x — iJ y + ifij — fi 2 . (2) 

H + ji 2 is positive definite, and if /i = m, with m S 
{-J, -(J — 1), . . . , (J — 1), J}, J = N/2 being the total 
angular momentum, it possesses a non degenerate ground 
state with E — obeying 

(aJ x - iJ y + imj |*) =0. 

Since this equation is linear the ground state can be easily 
constructed, which yields 

|*) = A/"(7, m) exp[-7 J z ] \m y = m) (3) 

where we have introduced the real anisotropy parameter 
7 through tanh(7) = a > 0. It is interesting to note that 
an anisotropy in the spin coupling is reflected here in the 
non-unitary term exp{— "fJ z }- In the fully isotropic limit 
7 = 0, the ground state is the state m y = m) which is 
entangled for all \m\ < J. In the maximally anisotropic 
case 7 = oo, the ground state is \m z = — J), which is a 
product state. The loss of entanglement in this case is 
due to the non- unitary term exp{— 7</ z }. Changing 7, 
e.g. as function of time, from 00 to causes a smooth 
transition from a factorized to an entangled many-body 
state. 

Due to the symmetry of the coupling all matrix ele- 
ments of the Hamiltonian between states corresponding 
to different total spin J vanish exactly even for time- 
dependent parameter. Thus even though the ground 
state of ^ becomes degenerate for 7 = with respect 



to the total spin J [13, the system cannot undergo a 
first order quantum phase transition upon changing 7. 
In fact the degeneracy in J at 7 = can easily be lifted 
by adding a term — AJ 2 to (Q, which has no effect on 

In the following we will restrict ourselves to the most 
interesting special case m — 0. As has been shown in 
[2l| and [23 the collective state \m = 0) has the largest 
global entanglement and should thus be considered as 
the state with maximum V-particle entanglement. A 
generalization to arbitrary m values is rather straight 
forward but less instructive. An additional feature of 
the m = case is the presence of a super-symmetry of 
the LMG Hamiltonian O. As a consequence m every 
spin sector J the spectrum of Q has for all values of 
7 a nondegenerate ground state and all excited states 
are pairwise degenerate [l5l | . As shown in [l6| the energy 
gap between the ground state and the pair of first excited 
states does not close. 

For m — the ground state 10 reads explicitly 

with Pj being Legendre polynomials. 

III. BI-PARTITE ENTANGLEMENT 

In the following section we discuss the entanglement 
between two arbitrary partitions of the N particle sys- 
tem in the SUSY ground state of the LMG model. As 
mentioned in the introduction it is not important here 
how the partitioning is done. Due to the symmetry of the 
Hamiltonian only the number of particles in each parti- 
tion is of relevance. 



A. Entropy of entanglement and distribution of 
Schmidt coefficients 

A generally accepted quantitative measure for the en- 
tanglement between two subsystems 1 and 2 if the total 
system is in a pure state | "J) is the von Neumann entropy 
of either of the two subsystems (entropy of entanglement) 

= -tri{pi lnpi} = -tr 2 {/92hi/92}, 

where 

P1.2 =tr 2 ,i{|tf)(tf|}, 

are the reduced density matrices. S(^) is essentially a 
measure for the information loss due to division of the 
system and ignoring one of the subsystems. If there is 
entanglement between 1 and 2 in the original pure state 
1 4") of the total system, the entropy is nonzero. On the 
other hand if |\E') factorizes there is no information loss 
if we ignore one subsystem and the entropy vanishes. 
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The von Neumann entropy for pure states S(JS) is iden- 
tical to the minimum relative entropy of entanglement 
-E^O^O H^with respect to all bi-partite separable states 
a e S 2 



Pi > o, 



i / 



Pi 



i. 



E 2 (^) = minSm\a), 
S(y\\a) = tr(plog 2 p- plog 2 CT) 



(5) 
(6) 



and p = |*)(4 f |. 

Calculating the von Neumann entropy of a many- 
particle system is in general a very nontrivial task due 
to the exponential growth of the relevant Hilbert space. 
We will show now that the von Neumann entropy can be 
related to the variance of the distribution of Schmidt co- 
efficients, arranged in an appropriate order, in the limit 
of a large number of particles. For the symmetric spin 
states considered here this variance can easily be calcu- 
lated, which will be done in the following subsection. 

Let | denote a pure state of a quantum system con- 
sisting of two parts labeled 1 and 2. In the case of finite 
dimensional spaces, Schmidt's theorem [24| asserts that 
any state | VP) in the Hilbert space Hi ® H 2 can be written 
in the form: 



unknown coefficient. In order to calculate this coefficient 
it is obviously not sufficient to use log 2 x as an estimate 
for S(tf). 

It is possible, however, to find a better estimate for 
S(^) in terms of the variance of the distribution of ap- 
propriately ordered Schmidt coefficients. To show this we 
first note that, as shown in [25|, the following inequality 
holds 

/(a;)| 2 log 2 \f(x)\ 2 dx< i(l + log 2 ^e) + log 2 Ax 

(10) 

where Ax 2 = J(x — x) 2 \ f(x)\ dx is the variance of a 
probability distribution f(x). This inequality becomes 
an equality if f(x) is a Gaussian function. 

For large values of x the sum in eq.© can be written 
as an integral with A 2 ^ — > A 2 (to), representing a con- 
tinuous, smooth probability distribution. The functional 
form of this distribution depends on the ordering of the 
Schmidt coefficients A m . For symmetric states this order- 
ing can be chosen in such a way, that A 2 (m) becomes to 
a good approximation a Gaussian function |26|. In this 
case the entropy of entanglement is given by 

5'(*) = i(l + log 2 7re)+log 2 AA ) (11) 
where AA is the variance of the Schmidt coefficients. 



m— 1 



<i>< 



$(2) 



(7) 



where x ^ min{di,d 2 }, d\ and d 2 being the dimen- 
sions of the corresponding Hilbert spaces, j and 

| $[ l 2 ' ) ^| are sets of orthonormal states for the sub- 
systems 1 and 2, respectively, and A m are the positive 
Schmidt coefficients obeying the sum rule 



E A * 



i. 



(8) 



It is easy to see that the entropy of entanglement is re- 
lated to the Schmidt coefficients via 



S(V)=-J2 A m log 2 X 2 



(9) 



The Schmidt rank, i.e. the number of Schmidt coeffi- 
cients x provides a simple upper bound for the entropy 
of entanglement S(^) < log 2 x- If X ~ min{di, d 2 }, i.e. 
if it scales exponentially with the number of particles, 
log 2 x is a polynomial function of the number of particles 
in the smaller of the two subsystems. For the symmetric 
coupling considered here, the dimension of the relevant 
Hilbert space d increases only linear in the number of 
particles implying a logarithmic scaling of log 2 x with 
the system size. Thus also the von Neumann entropy 5* 
is expected to scale logarithmically, however with a yet 



B. Clebsch-Gordan decomposition and bi-partite 
entanglement 

We will now calculate the entropy of entanglement us- 
ing ea. l|ll[l by finding a suitable bipartite decomposition 
of the ground state The scaling of S(^>) with the 
system size will be studied in detail and in particular the 
prefactor of the logarithm determined. In the limit of a 
totally isotropic spin coupling 7 = an explicit analytic 
expression for the variance of the Schmidt coefficients 
and the entanglement can be given. For nonzero values 
of 7 numerical results will be presented. 

We start from equation (0J and make use of the 
Clebsch-Gordan decomposition of total angular momen- 
tum J = N/2 into J 1 and J 2 with J\ + J 2 = J: 

|*)=^( 7 )e-^K=0) (12) 
= Mb) E E C ^,ti <oH7) \Ji mi) ® I J 2 ma) 

Here C^^^ are the Clebsch-Gordan (Wigner) coeffi- 
cients and to = toi + to 2 . d J t (/3) are the rotation 
matrices defined as H3 



d J m , rn ((3) = {Jm'\e-^ J «\Jm) 



(J + m)UJ-m)\ ( . (3 

-T TT7 -77 Sin — 

(J + m')!(J - m')! V 2 



cos — 
2 



(13) 

m-\-m ! 



X P 



(m—m , m+TTi ) 
J—m 



(cos/3), 
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where P^^^x) are Jacobi polynomials. Since we are 
considering the special decomposition for J = Ji + J 2 , 
the Wigner coefficients have the binomial distribution: 



tflJi J 2 J\ 
\ mim^m / 



2J X 
Ji + mi 



2J 2 

J2 + f»2 



2J 

J + to 



(14) 



This relation combined with l|13|) gives the following de- 
composition for the wave function |\&): 



Ji J2 

\*)= J2 V, m2 (7)|Jimi) 

mi=— Ji ms=- J2 

with 

A mi , m2 ( 7 )=^(7) (0 mi+m2 x 



|J 2 m 2 )(15) 



(J!) 2 (2J!)!(2J 2 )! 



(2J)!( Jx + m 1 )!(J 1 - mi)!(J 2 + m 2 )!(J 2 - m 2 )\ 



(oothl) 



j(— mi— m.2, mi+m2) 
J 



(cosh 7). (16) 



(|15JI is separable if all coefficients A mitrrl2 factorize. This 
is the case in the limit 7 — * 00, where the coth term in 
(|16[) approaches unity and the Jacobi polynomials factor- 
ize in mi and to 2 . 



different values of the total spin J = N/2. For J\ <C J a 
logarithmic scaling with prefactor 1/2 is evident. When 
Ji approaches J/2 the entropy saturates since S is sym- 
metric with respect to the replacement J\ <-> J — Ji . It 
is important to note that for J\ <C J the von Neumann 
entropy S does not depend on J. 




FIG. 1: Entropy of entanglement for 7 = as function of 
the logarithm of the subsystem spin Ji for different values of 
the total spin J. Due to the symmetry of S with respect to 
Ji <-> J — Ji the curves saturate at Ji = J/ 2. 



isotropic spin co 



2. anisotropic spin coupling 7 / 



Making use of the results of subsection A the summa- 
tion in (|12f> can be carried out explicitly for the limit of 
isotropic spin coupling. In this case one sees from l|16|) 
that only coefficients A mi ^ m2 with mi + m 2 = survive. 
Thus one has the following decomposition 



l*)= E G ii 



j 2 J1+J2 

m 



I Ji m) <g> I J 2 - m) 



m=— Ji 



where J 2 > Ji was assumed without loss of general- 
ity. The Schmidt coefficients are, therefore, the Clebsch- 
Gordan coefficients. In the limit J 2 3> Ji , they have a 
Gaussian form |26j 



Jl J 2 J1+./2 



exp 



where |to| < Ji <C J 2 . The Gaussian form of the coeffi- 
cients allows to make use of relation (|ll|l to calculate the 
von Neumann entropy for J\ <C J: 

S~ ilog 2 J 1 + i(l + log 2 7re), 7 = (17) 

In fig^ we have plotted the von Neumann entropy for 
7 = as a function of the subsystem spin Ji = iVi /2 for 



If the spin coupling is anisotropic, i.e. if 7 ^ 0, the 
double sum in cqs. (|12fl or (|15f) remains. Thus in order 
to discuss the influence of a finite 7 it is necessary to 
explicitly evaluate the sum in @. We have done this 
numerically for a total particle number up to 200 and 
subsystems up to 100 particles. The results are shown in 
figs|21and|3] As can be seen from fig[21in contrast to the 
isotropic case 7 = 0, the entropy is no longer independent 
on the total spin if 



7 > Tcrit 



1 
J 



In the thermodynamic limit the critical point is 7 = 0. 

As can be seen from fig |3] for any 7 > J -1 the entropy 
becomes a function of the logarithm of the fraction of 
particles Ji/J = A^/AT. 

Our numerical calculations suggest for J\ <^ J and 
7 3> J~ 2 a functional dependence of the form 

5- /(7)log 2 ( Ji/j), 7>7crit 
The reduction of entanglement with increasing 7 is ex- 



pected. The state A/"(7)e 



0) is maximally en- 



tangled for 7 = and the prefactor e ' yJz corresponds to 
a local non-unitary operation which always decreases the 
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FIG. 2: Entropy of entanglement for J = 100 (full line) and 
J = 200 (dashed line) as function of logarithm of subsystem 
spin Ji for different values of 7. One recognizes that in con- 
trast to the isotropic case the entropy now depends on the 
total spin J. 

5, 1 1 1 < 1 1 1 
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-8 -7 -6 -5 -4 -3 -2 -1 

log 2 J/J 

FIG. 3: Entropy of entanglement for different values of 7 
as function of log 2 (Ji/J) for J = 100 (solid line) and J = 
200 (dashed line). For 7 J > 1 the curves become virtually 
indistinguishable . 
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FIG. 4: Ordered distribution of normalized Schmidt numbers 
for different values of 7J. 



There is no obvious distinction of the point 7 = J -1 
in the properties of the system. The system does not un- 
dergo a phase transition at this point. Due to the SUSY 
the qualitative structure of the spectrum is the same for 
all values of 7 and there is alwas an energy gap between 
the ground and first excited state. Thus the question re- 
mains whether there are any physical signatures for the 
change of the scaling behavior of the entanglement at 

7 = 7crit- 



IV. GEOMETRIC ESTIMATE FOR GLOBAL 
ENTANGLEMENT 

In the previous section we have discussed the bi-partite 
entanglement of two partitions of the N spin 1/2 system. 
We have seen (cf. fi^3Jl that for 7 > 7 cr j t the von Neu- 
mann entropy has a maximum value independent on J. 
In this section we will quantitatively analyze this max- 
imum by determining the iV-partite or global entangle- 
ment Em of the SUSY ground state (@J which is an upper 
bound to the bipartite entanglement i?2- Although it is 
not possible to obtain an analytic expression for En, we 
can determine a very good estimate for it given by the 
geometric measure of entanglement. 



amount of entanglement. For large values of 7 the state 
becomes eventually separable. 

The most peculiar feature of the von Neumann en- 
tropy is the change of the scaling behavior with J\ from 
S ~ log 2 Ji for 7 < 7 crit to S ~ log 2 Ji/J for 7 > 7 crit . 
The role of 7 cr it = J -1 and the change of the scaling 
behavior is also reflected in the distribution of ordered 
Schmidt numbers. As can be seen from figQJthe fall-off 
of the Schmidt numbers A m becomes exponential when 
7 J exceeds unity. 



A. Relative entropy and geometric measure of 
entanglement 

A many-particle state is called iV-partite separable, if 
it can be written as a product of states of all N sub- 
systems. Obviously a bi-partite entangled state is al- 
ways iV-partite entangled but not vice versa. A quanti- 
tative measure of many-particle or global entanglement 
of a state p is the minimum relative entropy which deter- 
mines the minimum distance between p and the set Sn 
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of TV-partite product states a [2j| 

E N = min S{p\\a) (18) 

ct£5jv 

where 

S{p\\cr) =tr(/3log 2 /3-/3log 2 cr), (19) 
cr G 5/v being an Y-partite separable state 

N 

0- = 5^PiPl® l4® -<8> Pn ( 2 0) 
z=l 

with ^ > 0, and JYPi — !• For the bi-partite case £W is 
equivalent to the entanglement of formation |23j |. which 
in the case of pure states is identical to the von Neumann 
entropy. 

Since the set Sn is smaller than iS 2 for any partioning, 
Sn C ^2, it follows immediately 

E N {V) > E 2 (V), 

i.e. the global entanglement represents an upper bound 
to the bi-partite entanglement. 

In order to compute En for any state p, one has to 
find its closest product state a. This is in general a quite 
difficult task and can be done only in very special cases. 
There is however a lower bound to En which gives a 
good estimate for the behavior of the global entangle- 
ment. This lower bound is the geometric entanglement 

mm. 



where 



E G {9) = -21og 2 A max (4-), 



A max (*) = max I I 

<t> 



(21) 



(22) 



is the maximum overlap of I*) with an TV-partite sepa- 
rable state \<fi). E G (^) is not an entanglement monotone 
and thus in the strict sense not a valid measure of entan- 
glement. It does give however a close lower bound to En 
which for some states such as the Dicke states is a tight 
bound, i.e. E G = En HJ- The geometric entanglement 
can easily be calculated for states which are permutation 
invariant, which is the case for the SUSY ground state 

Q. 



B. Geometric measure of entanglement for the 
SUSY state 



Calculating the overlap of \<fi) with Q and maximizing 
it with respect to the real parameters a, (3 and £ leads to 

10) = \m z = -J). 

The corresponding entanglement eigenvalue reads 



A max (7) 



2 J J! y/Pj (cosh 2 7 )' 



(23) 



1. isotropic spin coupling 7 = 
In the isotropic case Il2:il) reduces to [2l| 

t (7 = 0) 



A,: 



2 J J\ 



(24) 



and thus the geometric entanglement is given by 
E G {^) = \ log 2 J. Since the SUSY state for 7 = is 
the Dicke state | J, m y — 0) the geometric entanglement 
is a tight lower bound to the relative entropy and thus 



B Ar (*) = B G (*) = iiog 2 J. 



2. anisotropic spin coupling 7 7^ 



(25) 



It is obvious that the largest entanglement is obtained 
for 7 = where the maximum overlap A max with separa- 
ble states is smallest. On the other hand for 7 — > 00 the 
state becomes identical to the separable state \m z = — J). 
The same conclusion can of course be obtained from 
eq. (|23|l employing the asymptotic expansion of the Leg- 
endre polynomials. 
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To calculate the geometric entanglement Eq or equiv- 
alently the maximum overlap A max of the SUSY ground 
state (0} with A^-partite separable states it is sufficient 
to construct the most general A^-partite separable state 
which is invariant under permutation of spins 21]. This 
state is given by rotations of the state \m z = — J): 

|0) = e -fai. e -^i» e -*€i. \ mz = 



FIG. 5: Geometric measure of entanglement as function of 
log 2 J for different values of the anisotropy parameter 7. One 
recognizes a saturation at J > 7~ it . 

In fig|5]we have plotted the geometric entanglement as 
a function of J for different values of 7. For sufficiently 
small values of J one recognizes a logarithmic growth 
which saturates when J exceeds the value 7 -1 . One can 
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easily obtain an analytic expression for the saturation 
value of Eq. Making use of the asymptotics of the Leg- 
endre polynomials for large J and 7 7^ 



Pj (cosh 27) 



1 



(2J-1) 



largo J VI -e-^ " 2 J J! 
one arrives at the simple expression 



A max ( 7 ) - (1 - e-^Y' A 



leading to 



E G (l,J) £^7 --log 2 (l- e - 4 ^). 

Comparing the numerical values for E2 obtained in the 
previous section one finds that Eq < E 2 . This shows that 
in the case of non-isotropic coupling 7 7^ the geometric 
entanglement is not a tight lower bound to the global 
entanglement, i.e. here Eq < Ejf. Thus Eq can only be 
used as qualitative measure for the global entanglement. 



V. CONCLUSIONS 

In the present paper we have studied the bi-partite 
entanglement between blocks of spins in the anti- 
ferromagnetic Lipkin-Meshkov-Glick model, which de- 
scribes a nonlinear coupling of collective spins, under con- 
ditions of super-symmetry. The super-symmetry of the 
model allows for an explicit construction of the ground 
state which undergoes a smooth transition from a sep- 
arable to a maximally entangled state when changing 
the anisotropy of the collective spin coupling. Making 
use of the Clebsch-Gordan decomposition of angular mo- 
menta, the von Neumann entropy which quantifies the 



bi-partite entanglement can be calculated analytically in 
the isotropic case or numerically in the case of anisotropic 
coupling. Although the structure of the spectrum stays 
always the same with one nondegenerate ground state 
and pairwise degenerate excited states, and no level cross- 
ing or merging occurs, the entanglement shows a dis- 
continuous behavior at the isotropy point. When the 
anisotropy parameter 7 vanishes exactly, the von Neu- 
mann entropy grows logarithmically with the number of 
particles in the subsystem. For any nonvanishing value of 
7 (in the thermodynamic limit) the entropy saturates at 
a finite value determined by 7. The maximum bi-partite 
entanglement can be estimated by the geometric measure 
of global entanglement, which has been determined ana- 
lytically. Furthermore in this case the entropy becomes 
a function of the ratio of particle number in the subsys- 
tem to the total particle number rather than a function 
of the subsystem size alone. For finite systems the tran- 
sition between the two cases happens at a small but fi- 
nite value of 7 corresponding to the inverse of the total 
number of spins. A discontinuous scaling behavior of en- 
tanglement is usually attributed to level crossings and 
quantum phase transitions. This is not the case in the 
present system and the question remains whether there 
are any physical signatures of the discontinuous transi- 
tion. 
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